Chiral and deconfinement transition from Dyson-Schwinger equations 
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We determine the quark condensate and the dressed Polyakov loop from the finite temperature 
Landau gauge quark propagator evaluated with f/(l)-valued boundary conditions in an approxima- 
tion to quenched QCD. These gauge invariant quantities allow for an investigation of the chiral and 
deconfinement transition. We compare results from Dyson-Schwinger equations on a lattice with 
infinite volume continuum results and study the temperature and quark mass dependence of both 
quantities. In particular we investigate the chiral condensate and the dressed Polyakov loop in the 
chiral limit. We also consider an alternative order parameter for the deconfinement transition, the 
dual scalar quark dressing, and compare with the dressed Polyakov loop. As a result we find only 
slightly different transition temperatures for the chiral and the deconfinement transition at finite 
quark masses; in the chiral limit both transitions coincide. 

PACS numbers: 12.38.Aw, 12.38.Lg,11.10.Wx 



I. INTRODUCTION 

The chiral and the deconfinement transition of QCD 
are a subject of continuous interest both, from a theoreti- 
cal and from an experimental point of view. Ongoing ex- 
periments at the Relativistic Heavy Ion Collider (RHIC) 
and future runs at LHC and FAIR probe the transition 
between the high and low temperature states of strongly 
interacting quark matter. The transition is characterized 
by the breaking and restoration of chiral and center sym- 
metry. Hereby it is well known that both transitions are 
phase transitions in the strict sense of the word only in 
opposite limits of the theory. In the chiral limit m = 
chiral quantities such as the chiral condensate (ipip) or 
the scalar part -B(O) of the inverse quark propagator at 
zero momentum are order parameters of the chiral phase 
transition. In the absence of dynamical quarks, m — > oo, 
pure gluonic systems are center symmetric in the low 
temperature phase. The deconfinement phase transition 
is then related to center symmetry breaking indicated by 
a non vanishing order parameter as e.g. the Polyakov 
loop L fl. In the real world, quarks carry (light) masses 
and both transitions are crossovers as nicely established 
by lattice QCD, see e.g. @,SH@ 

and references therein. 
The theoretical search for an underlying mechanism 
that links quark confinement and chiral symmetry break- 
ing has brought many interesting results in the past 
years. Topological field configurations such as center 
vortices have been argued to connect both phenomena, 
see e.g. for a review. In 0] a close connection be- 
tween center vortices and the spectrum of the Dirac oper- 
ator has been revealed; removing those vortices from the 
gauge field ensemble resulted in a vanishing string ten- 
sion and a zero chiral condensate via the Casher-Banks 
relation. 

Very recently a number of works @, i, [H OH, 0, GJ, 
[Til . [l5| have explored the relation between chiral and 
deconfinement signatures in the spectrum of the Dirac- 



operator. In [ll| it has been shown that the low ly- 
ing eigenmodes of the Dirac operator are responsible for 
both, chiral symmetry breaking and confinement: spec- 
tral sums of the Dirac operator may be truncated to 
the lowest lying modes and still lead to a linear rising 
quark-antiquark potential extracted from the Polyakov- 
loop correlator. 

One of these spectral sums, the dual condensate or 
dressed Polyakov loop, is of particular interest since it 
relates the chiral condensate to the Polyakov loop. This 
quantity has been studied in lattice gauge theory [l2| 
and with Dyson-Schwinger equations on a torus [16[ with 
qualitatively similar results. Both of these studies dealt 
with the dressed Polyakov loop on a compact manifold 
at finite quark masses. In this work we add more details 
to these investigations, in particular considering the infi- 
nite volume limit, the chiral limit and the large tempera- 
ture behavior of these quantities. We also investigate an- 
other order parameter for the deconfinement transition, 
the dual scalar quark dressing function. Compared to 
Ref. [T(| we also modified the renormalization condition 
of the quark mass function in the torus formulation, thus 
eliminating temperature effects in the bare quark mass. 
As a result we now obtain almost similar transition tem- 
peratures for the chiral and the deconfinement transition 
for quark masses roughly corresponding to an up-quark 
and identical transition temperatures in the chiral limit. 

This work is organized as follows: In section |TT] we 
outline our two order parameters for the confmcmcnt- 
dcconfmcmcnt transition. Here we consider (i) the dual 
quark condensate as defined in Ref. [l2j and (ii) the dual 
scalar quark dressing function as defined below. Both 
quantities share the property that they can be extracted 
from the momentum dependence of the quark propagator 
and are therefore accessible with functional methods. In 
this work we use the Dyson-Schwinger equations (DSEs) 
for the quark propagator in Landau gauge to determine 
these order parameters. The temperature dependence 
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of the propagator is evaluated in the Matsubara formal- 
ism. Wc employ a truncation scheme for the quark-DSE 
that uses 5f/(2)-lattice results for the temperature de- 
pendent gluon propagator [T3| and a temperature depen- 
dent ansatz for the quark-gluon vertex as input. This 
scheme is detailed in section IIII1 where we also discuss 
our numerical procedure. In section IIVI we present our 
numerical results in the infinite volume/continuum for- 
mulation of the DSEs and compare with the results of 
Ref. [l6[ obtained from a formulation of the DSEs on 
a torus. Our infinite volume approach also enables us 
to explore the chiral limit of the theory. We compare 
results for the order parameters specified in section |TT] 
and comment on relations to the chiral phase transition. 
Our results are summarized in the concluding section Ivl 
Throughout the paper we work in SU (2)- Yang-Mills the- 
ory, leaving the SU(3) case for future work. 



II. ORDER PARAMETERS FOR 
DECONFINEMENT 

Quantities that transform non trivial under center 
symmetry and thus qualify as order parameters for de- 
confinement are in general not easily accessible by func- 
tional methods. However, there has been considerable 
progress on this issue in the past years. In Ref. [HI, the 
ghost and gluon propagators of Landau gauge Yang-Mills 
theory have been used to determine the Polyakov loop po- 
tential within an effective action approach. In this work 
we will focus on an alternative strategy that allows to 
extract chiral and deconfinemcnt order parameters from 
a (generalized) quark propagator. This approach estab- 
lishes links between dynamical chiral s ym metry breaking 



a sum over closed loops winding n-times around the com- 
pact time direction: 



and quark confinement in 



ncal chiral s ym metry break 
QCD. @, i, Mi, 111, [H. 



A. The dual condensate 

The dual quark condensate as order parameter for cen- 
ter symmetry has emerged from a series of works connect- 
ing spectral sums of the Dirac propagator with Polyakov 
loops and their correlators [1, 0, [l(| ■ Within the frame- 
work of lattice gauge theory the dual condensate S n has 
been introduced in Ref. [l2|]. It is defined by the phasc- 
Fouricr-transform 
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of the ordinary quark condensate evaluated with respect 
to a phase e l</3 in the temporal direction of the Euclidean 
theory. This phase is introduced by the generalized, 
£/(l)-valued boundary condition tp(x, 1/T) = e^tp^.O) 
in the temporal direction. For the usual antiperiodic 
boundary conditions for fcrmions we have tp = tt, whereas 
ip = corresponds to periodic boundary conditions. The 
(/^-dependent quark condensate (ipip)<p ' 1S proportional to 
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Here C denotes the set of all closed loops I with length 
|/| on a lattice with lattice spacing a. Furthermore m is 
the quark mass. U (I) stands for the chain of link vari- 
ables in a loop I multiplied with appropriate sign and 
normalization factors, see Ref. [H] for details. Each loop 
that closes around the temporal boundary picks up fac- 
tors of e ±Iip according to its winding number n(l). The 
Fourier transform in (JTJ) projects out those loops with 
winding number n. The dual condensate T,i then corre- 
sponds to loops that wind exactly once and is called the 
'dressed Polyakov loop' This quantity transforms 

under center transformation identically as the conven- 
tional Polyakov loop [l[ and is therefore an order param- 
eter for center symmetry. The numerical agreement be- 
tween dressed and conventional Polyakov loop has been 
established for gauge groups SU(3) [13| and, remarkably, 
also for the center less G(2) [3]. 

These lattice calculations have been carried out in 
quenched QCD and in a formulation where the nontrivial 
[/(l)-valued boundary conditions are introduced for the 
measured operators only. Such a formulation avoids the 
complications associated with the introduction of non- 
trivial boundary conditions into the generating functional 
of QCD with fermions. The latter option is formally 
equivalent to the introduction of an imaginary chemical 
potential leading to the pattern of Roberge- Weiss peri- 
odicity of the generating functional [2(| . In this work we 
will approximate the quenched lattice formulation and 
therefore do not encounter the Roberge- Weiss symmetry. 
A different approach is followed, however, in the (un- 
quenched) renormalization group framework of Ref. [l9( , 
where the Roberge- Weiss symmetry is explicitly taken 
into account. 

Furthermore, a comment on regularization is in or- 
der here. Below wc will determine the dual conden- 
sate Ei from Dyson-Schwingcr equations in the infinite 
volume and continuum limit. For non- vanishing quark 
masses the quark condensate is quadratically diver gent 
and needs to be properly regularized, see e.g. Ref. [21| 
for a recent discussion. Correspondingly, the loop expan- 
sion Eq. ^ breaks down for a — ► 0. In this work we will 
employ an ultraviolet regulator in the form of a simple 
cutoff. For large enough cutoffs the regulator dependent 
part of the condensate is independent of temperature and 
therefore does not affect the chiral transition tempera- 
ture. The regulator dependent part of the condensate is 
also independent of the boundary angle. Consequently 
it does not appear in the dual condensate and does not 
affect the deconfinement transition temperature Td ec ei- 
ther. The latter property can also be demonstrated a 
posteriori by comparing T^ec from the dual condensate 
with the one obtained from the dual scalar quark dress- 
ing function introduced in the next subsection. Since 
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FIG. 1: The Dyson-Schwinger equation for the quark propa- 
gator. Filled circles denote dressed propagators whereas the 
empty circle stands for the dressed quark-gluon vertex. 



is an order parameter for the confinemcnt-dcconfincmcnt 
transition. This can be shown in a similar fash- 
ion as for the spectral sums discussed in [ll[: Using 
B(0,ui p (0,ip)) = l/4tr[5- 1 (0,w p (0,v3))] we have 



-up 



l/T 

d 3 x { dxotr&xolD-^O)- 1 (5) 



the (dual) scalar quark dressing function does not suf- 
fer from regularization problems, agreement between the 
corresponding transition temperatures indicates that our 
regularization procedure is adequate. This is indeed the 
case as demonstrated in section IIVCI 



B. The dual scalar quark dressing function 

Another order parameter for the deconfincment tran- 
sition is the phase-Fourier-transform of the scalar quark 
dressing function evaluated at lowest Matsubara fre- 
quency and zero momentum. To derive this quantity note 
first that in momentum space the Z7(l)-valued boundary 
conditions introduced above result in Matsubara modes 
ujp(nt, tp) = (2irT)(nt + p/2tt) in the ^-direction, which 
depend on the boundary angle ip € [0, 2tt[. The inverse 
quark propagator can then be written as 

S~ l {p 1 uj p ) = ij 4 u p C(p,u> p ) +ijiPiA(p,ujp) +B(p,Lj p ) , 

(3) 

with vector and scalar quark dressing functions C, A, B. 
A further tensor component proportional to a^ v is pos- 
sible in principle but can be omitted in all practical cal- 
culations p2| . 

The scalar quark dressing function B{p,u) v ) evaluated 
at p = and ui p (0, tt) = ttT is an order parameter for the 
chiral phase transition in the chiral limit. It also serves 
as indicator for the chiral crossover at non- vanishing bare 
quark masses in a similar way as the chiral condensate 
does. Its dual 



where D v is the (massive or chiral) Dirac operator evalu- 
ated under presence of the [/(l)-valued boundary condi- 
tions. A center transformation on introduces an ad- 
ditional phase factor z = e j27rfc /JV w j^ u k = 0, . . . , N c — 1 
for N c colors, which adds to the phase e lip for our U(l)- 
valucd boundary conditions. We then obtain 



2tt l/T 
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l/T 



^e-* B(O,w p (O,¥0) 



(4) 



/ ^ e ~ iV J d * x I dxotr&XolD-^O)- 1 ,(6) 





i.e. the dual scalar quark dressing Eg transforms un- 
der center transformations exactly like the conventional 
Polyakov loop and therefore acts as order parameter for 
the deconfinement transition in the heavy quark limit. 
This will be confirmed by our numerical results given be- 
low. 



III. DYSON-SCHWINGER EQUATIONS FOR 
THE QUARK PROPAGATOR 

A. The truncation scheme 

The Dyson-Schwinger equation for the quark propaga- 
tor Eq. © is displayed diagrammatically in Fig. [TJ At 
finite temperature T it is given by 



S 1 {p,uj p ) = Z 2 S Q 1 (p,uj p ) - C} 



Z 2 Z\ 



d 3 k 

.—^t S(k,uj k )T u (k,uj k ,p,ujp) Dp V {p- k,uj p - w k ) 
( 27T ) 



(7) 



r 



Here Dp U denotes the (transverse) gluon propagator in 
Landau gauge and we have introduced a reduced quark- 



gluon vertex T v , by defining T^jf 1 = ig^-Y v . 
quark propagator is given by S^ 1 {p) = vy -p 4 
where m(/i 2 ) is the renormalized current quark mass 



The bare 
Z m m(jj, 2 ), 



The Casimir factor Cp = (N 2 — l)/N c stems from the 
color trace; in this work we only consider the gauge group 
SU(2). The wave function and quark mass renormaliza- 
tion factors, Z 2 and Z m , are determined in the rcnor- 
malization process. The ghost renormalization factor Z% 
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FIG. 2: Quenched 5(7(2) lattice results for the transverse dressing function Zr(q) and the longitudinal dressing function Zl(q) 
of the gluon propagator [l7| together with the fit functions Eq. (I10|l . Note that the ultraviolet part of the lattice data has not 
been corrected for hyper-cubic artefacts; the fits (straight lines; see text) are adapted to the corresponding corrected data [25(] 
at temperature T = 0. 



is canceled by a corresponding factor in our model for 
the quark-gluon vertex discussed below. Furthermore we 
used Z\ — 1 for the renormalization factor of the Landau 
gauge ghost-gluon vertex. The quark dressing functions 
A(p,ujp), B(p,u>p) and C{p,uj p ) can be extracted from 
Eq. ||7J) by suitable projections in Dirac-spacc. 

In order to solve this equation we have to specify ex- 
plicit expressions for the gluon propagator and the quark- 
gluon vertex. At finite temperatures the tensor structure 
of the Landau gauge gluon propagator contains two parts, 
one transversal and one longitudinal to the heat bath. 
The propagator is then given by (q = (q,u) q )) 

DM = Z -^pIM + ^p l ,M (8) 

with transverse and longitudinal projectors 

p^(q) = *V(?)--P£(?), (9) 

with (i,j = l... 3). The transverse dressing ZT(q,w q ) is 
also known as magnetic dressing function of the gluon, 
whereas the longitudinal component Zl (q) is called elec- 
tric dressing function of the gluon propagator. At zero 
temperatures Euclidean 0(4)-invariance requires both 
dressing functions to agree, i.e. Zr(q) = Zt{q) = Z(q). 

For the momentum range relevant for Eq. (|7|) we nowa- 
days have very accurate solutions for the gluon dress- 
ing function Z(q) at zero temperature from both, lattice 
calculations and functional methods, see e.g. [23| and 
references therein. The temperature dependence of the 
gluon propagator, however, is much less explored, see 
e -g- [HI- In [13 a combined study on the lattice and 



from Dyson-Schwingcr equations records a very differ- 
ent temperature dependence of the electric and magnetic 
parts. This is shown in Fig. [2J The temperature effects 
on both the magnetic and electric dressing functions are 
such that there are almost no effects when comparing the 
T = result with T = 119 MeV. Further increasing the 
temperature to T = 298 MeV and T = 597 MeV sig- 
nificantly decreases the bump in the magnetic dressing 
function around p 2 = 1 GeV 2 . There is no indication 
that this decrease takes special notice of the critical tem- 
perature T c rts 300 MeV for quenched QCD with gauge 
group SU (2). The opposite seems to be true for the elec- 
tric part of the propagator. Here from T = 119 MeV to 
T = 300 MeV one observes a clear increase of the bump 
in the dressing function Zl (q) and a subsequent decrease 
when the temperature is further raised to T — 597 MeV. 
Pending further investigation it seems reasonable to as- 
sume that the maximum of the bump is reached at or 
around the critical temperature T c sa 300 MeV. 

Although the lattice data still have considerable sys- 
tematic errors (l~7| they may very well correctly repre- 
sent the qualitative temperature dependence of the gluon 
propagator. We therefore use a temperature dependent 
(qualitative) fit to the data as input into the DSE; this fit 
is also displayed in Fig. [5] (straight lines). Note that for 
momenta around p 2 = 6 GeV 2 the lattice data show con- 
siderable systematic errors due to hyper-cubic artefacts. 



T[MeV] 





119 


298 


597 


a T (T) 


1 


1 


1.34 


1.65 


a L (L) 


1 


1 


0.8 


4.0 



TABLE I: Temperature dependent fit parameter of Eq . (|10fl . 
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For the fits we have eliminated these artefacts by com- 
parison with the T = lattice result of Ref. [25[ , where 
these errors have been corrected. The fit functions read 



(q,uj q ,T) 



2 A 2 



q 2 A 



('l 2 



A 2 



A 2 ) 2 |V<? 2 + A2a T ,L(T) 
A)«(m) ln[g 2 /A 2 + l] ^ 

47T 



(10) 



with the temperature independent scale A = 1.4 GeV and 
the coefficient c = 9.8 GeV 2 . For gauge group SU{2) we 
have Po = 22/3 and 7 = —13/22 in the quenched theory 
and we renormalize at a(fi) = 0.3. The temperature de- 



pendent scale modification parameters <zt,l (T) are given 
in table HI In order to extend this fit to temperatures not 
given in the table we assume clt,l(T) to be temperature 
independent below T = 119 MeV and only slowly rising 
above T = 597 MeV. For T G [119,597] MeV we use 
cubic splines to interpolate smoothly between the val- 
ues given in table HI We expect the systematic error of 
this procedure to be of the same order as the systematic 
errors inherent in the lattice data. We also inherit the 
scale determined on the lattice using the string tension 
0.44 GeV [13 ■ 

For the quark-gluon vertex with gluon momentum q = 
(q,u> q ) and the quark momenta p = (p,ui p ),k = (k,u>k) 
we employ the following temperature dependent model 



J 



T v (q, k,p) = Z 3 <5 4 „7 4 



C(k)+C(p) , s _ A(k)+A(p)\ ( di q 2 fp a(n)ln[q 2 /A 2 + l] 




r 



q 2 A 2 + q 2 



-in 



25\ 



(11) 



where 8 = —9/44 is the anomalous dimension of the 
vertex. Note that because of 7 + 25 = —1 the gluon 
dressing function together with the quark-gluon vertex 
behaves like the running coupling at large momenta; this 
is a necessary boundary condition for any model inter- 
action in the quark DSE. The dependence of the vertex 
on the quark dressing functions A and C is motivated 
by the Slavnov- Taylor identity for the vertex; it repre- 
sents the first term of a generalization of the Ball-Chiu 
vertex (2(| to finite temperatures. The remaining fit func- 
tion is purely phenomenological, see e.g. [27l ] where an 
elaborate version of such an ansatz has been used to de- 
scribe meson observables. Here we use di = 7.6 GeV 2 and 
c?2 = 0.5 GeV 2 . A moderate variation of these parame- 
ters shifts the critical temperatures of both, the chiral 
and the deconfincment transition but leaves all qualita- 
tive aspects of the results presented below unchanged. 



To our knowledge the truncation scheme defined above 
is the first of its kind within the DSE-approach that im- 
plements a realistic temperature dependence of the gluon 
propagator and the quark- gluo n vertex beyond simple 
ansaetze, see e.g. f22l [28l [29l |30| for previous approaches. 
The explicit expressions of the resulting DSEs for the 
quark dressing functions including Eq. (flTJ)) and Eq. pd|) 
are given in appendix [SJ Unless explicitly denoted oth- 
erwise all results presented in section ITVl of this work are 
obtained using these DSEs. 

For the sake of comparison, however, we also employed 
a simplified input, where all temperature dependencies 
are solely due to the quark propagator. To this end we 
fix the temperature dependent parameter a^t (T) in the 
gluon dressing function to their zero temperature values 
and remove the temperature effects in the vertex using 



J 



r„(q, k,p) = Z^ v 



di 



d 2 + q 2 A 2 



/3 a(/x)ln[g 2 /A 2 + l] 



4?r 



2S\ 



(12) 



r 



with (modified) fit parameters d\ = 10 GeV 2 and e?2 = 
0.2 GeV 2 such that the chiral transition temperature is 
the same as for the temperature dependent interaction. 

Finally, when investigating the effect of the details of 
the quark-gluon interaction on the phase transition an 
interesting question concerns the influence of the deep 
infrared behavior of the ghost and gluon propagators as 



well as the quark-gluon vertex on the phase transition. 
To this end we note that the expressions (TIT)]) and pd|) 
for the gluon propagator and the quark-gluon vertex em- 
ployed above correspond to decoupling type of infrared 
behavior in the sense specified in [23||. The correspond- 
ing expressions for scaling in the infrared are given in 
appendix [Cl where also corresponding numerical results 
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for the temperature dependent condensates are given. As 
a result we find that the chiral and deconfincmcnt tran- 
sitions is insensitive to the question of scaling vs. decou- 
pling in the deep infrared. 



B. Numerical procedure 

Infinite volume/continuum DSEs: 

The numerical solutions of the DSE are found using the 
conventional fixed point iteration method. The loop in- 
tegral and Matsubara sum is regularized using a sharp 
cutoff A = s/E x 10 2 GcV such that the integration and 
summation extends to momenta and frequencies with 
+q 1 < A 2 . This is requisite to restore 0(4) invariance 
at perturbative momenta as well as in the limit of zero 
temperature. We apply a MOM renormalization scheme 
which defines the renormalization constants in such a way 
as to fix the full propagator at a certain momentum. We 
explicitly verified that the resulting renormalized quark 
propagator is independent of the size of the ultraviolet 
cutoff by increasing/decreasing the cutoff by an order 
of magnitude. The renormalized quark mass m(/z) has 
contributions stemming from the explicit chiral symme- 
try breaking and from the quark condensate. Since the 
quark condensate is temperature dependent fixing the 
renormalized quark mass with temperature independent 
renormalization conditions means that the magnitude of 
explicit chiral symmetry breaking changes with temper- 
ature. In order to minimize this effect we choose a very 
large renormalization point such that the quark conden- 
sate term is irrelevant there and the temperature effects 
in m(/x) can safely be ignored. The renormalization con- 
ditions arc C(/i) = 1 and B(/i) = m(/i) with /i = {fl, irT) 
and /i 2 = 10545 GcV 2 . The choice m(/i) = 1.85 MeV 
then matches the renormalization conditions in Ref. [l6[ 
for mass 10 MeV at T = 200 McV and = 20 GeV 2 and 
corresponds roughly to an up quark mass. For \n t \ < 20 
the summation over Matsubara frequencies is performed 
explicitly The remaining sum is approximated by an in- 
tegral that can be done with standard Gauss-Legendre 
integration. 

In the chiral limit the transition temperature for the 
chiral phase transition is unambiguously defined by the 
temperature at which the symmetry is completely re- 
stored. This is signaled by the vanishing of the order pa- 
rameters i.e. the quark mass function at zero momentum 
or equivalently the vanishing of the chiral quark conden- 
sate. For explicitly broken chiral symmetry we encounter 
a crossover rather than a phase transition and the tran- 
sition temperature may depend on its definition. It can 
be defined by the peak of the chiral susceptibility 



X = d/(dm) (ipip) 
or the renormalized chiral susceptibility 
d 



(13) 



which can be made dimensionless by normalizing with 
T 4 for details see Refs. [HQ. An alternative definition of 
the transition temperature is via the (normalized) slope 
of the condensate with temperature 



sup 

T 



9 (jrtp) \ 
T 2 d T ) ' 



(15) 



The if dependent quark condensate is determined ac- 
cording to 



u) p (ip) 



(¥>)) (16) 



with the conventional quark condensate obtained for 
(p = 7r and multiplication with Z m . As mentioned in 
section Hi A[ the integral in (JT5J) is divergent for nonzero 
quark masses and needs to be regularized. To this end 
we introduce a cutoff at large enough momenta where 
the quark propagator in the integrand is temperature 
independent to very good approximation. The result- 
ing regularized condensate is not indicative as concerns 
its absolute value but retains its correct dependence on 
temperature T and boundary angle tp as long as the tem- 
perature is much smaller than the cutoff scale. This is 
the case for all results presented in the following. 

DSEs on a torus: 
In Ref. [l6[ first results for the ordinary and the dual chi- 
ral condensate have been presented in a formulation of 
the Dyson-Schwingcr equations on a torus. Such a formu- 
lation has the advantage of considerable simplifications 
as concerns the numerics but has the disadvantage of in- 
troducing artefacts due to the finite volume and lattice 
spacing. The calculations on a torus have been carried 
out at a three- volume of V = (5fm) 3 . This is sufficiently 
large to avoid significant volume effects even for quark 
masses as small as m(/2 2 = 20 GeV 2 ) = 10 McV used 
in Ref. For even larger volumes one obtains small 

changes of the order of a few percent in the chiral con- 
densate below the critical temperature and much smaller 
effects for larger temperatures; the value of the chiral 
transition temperature is unaffected. In addition there 
are effects due to the finite lattice spacing in temporal 
and spatial directions and the fact that the renormalized 
quark mass m = B(ttT, /2)/C(7rT, jl) has been treated 
as a constant with respect to variations in temperature 
[l6j ]. This is certainly correct for a very large renormal- 
ization point /I, but introduces temperature dependent 
artefacts for the value fx 2 = 20 GeV used in the torus cal- 
culation. In this work we remove these artefacts by the 
following procedure: we first solve the Dyson-Schwingcr 
equations in the infinite volume/continuum limit, i.e. in 
the formulation introduced in section IIIIl These equa- 
tions are renormalized at a very large renormalization 
point and consequently temperature effects in the renor- 
malization procedure are absent. From these results we 
read off the corresponding, temperature dependent renor- 
malization conditions at the small renormalization point 
jl 2 = 20 GcV 2 of our torus formulation. These are then 
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FIG. 3: Left diagram: Temperature dependence of the dressed Polyakov-loop £1 and the conventional quark condensate 
= {^V')^=7r together with their derivatives for m = lOMeV evaluated at two different ultraviolet cutoffs (thick lines and 
filled symbols for the large cutoff; thin lines and open symbols for the smaller cutoff). Right diagram: The same calculation in 
the infinite volume limit and at an extremely large ultraviolet cutoff. Note that the ordinary condensates in both plots have 
been shifted to agree at T = 200 MeV. 



applied to the calculations on the torus. This procedure 
is justified if and only if the truncation of the quark- 
DSE respects multiplicative rcnormalizability, which is 
the case for our truncation scheme. As we will see below, 
this change in the renormalization procedure leads to a 
significant change in the deconfinement transition tem- 
perature on the torus as compared to the one extracted 
in Ref. 0. 

The numerical details for the calculation on a torus 
are similar to the ones in the zero temperature limit, dis- 
cussed at length in Refs. 3l|. As concerns corrections for 



hyper-cubic artefacts we employ the 'cutting-the edges' 
prescription developed in [3l[ in the three spatial direc- 
tions of the torus. For the temporal direction we take 
into account as many Matsubara modes as necessary to 
arrive at approximately the same cutoff scale as for the 
spatial directions. This results e.g. in n t = 8 temporal 
Matsubara modes at T = 200 MeV. 

These prescriptions are adapted to minimize volume 
and discretization artefacts when comparing solutions 
from DSEs on a torus with the corresponding results in 
the infinite volume/continuum limit. At this point we do 
not claim that the remaining volume and spacing effects 
are similar to the corresponding ones in lattice calcu- 
lations. Nevertheless it is interesting to investigate the 
details of the volume and cutoff-dependence of our torus 
results systematically. While first results of this compar- 
ison are presented in the next section, a more detailed 
analysis is deferred to future work. 



IV. NUMERICAL RESULTS 

A. The infinite volume and continuum limit 

We first compare results from DSEs on a torus with 
results from the infinite volume/continuum DSEs at a 



quark mass that roughly corresponds to an up-quark. 
The resulting temperature dependent ordinary and dual 
condensates as well as their (normalized) temperature 
derivatives are shown in Fig. [3] The left hand diagram 
displays the torus results at two different ultraviolet cut- 
offs p m ax = 5.96, 7.45 GeV corresponding to lattice spac- 
ings of a = n/p max w 0.1, 0.08 fm, respectively. The right 
hand diagram shows our results in the continuum. Com- 
pared to Ref. [l6| we use a modified mass renormalization 
scheme on the torus as described in the last subsection. 
Furthermore we normalize the temperature derivatives 
with 1/T 2 to obtain dimensionlcss quantities. The last 
change affects both transition temperatures only by 1-2 
MeV. The effects of the modified mass renormalization 
procedure can be assessed by comparing our results with 
the ones shown in Fig. 3 of Ref. [lj| . We obtain a more 
pronounced chiral transition with a narrower peak in the 
temperature derivative of the ordinary chiral condensate. 
The corresponding transition temperature r w 300 MeV 
is unchanged. As for the dual condensate, or dressed 
Polyakov loop, the peak of the temperature derivative 
is shifted to smaller temperatures, i.e. from Td ec ~ 320 
MeV in Ref. to our value T dec w 308 MeV for both 
cutoffs. Despite this shift, the deconfinement transition 
remains separate from the chiral transition. 

This finding is also visible in the infinite vol- 
ume/continuum limit. Although there are still differences 
in the details of the temperature dependence of the ordi- 
nary and dual condensate between torus and continuum 
DSEs, the related peaks in the temperature derivatives 
are similar; here we find r « 301 MeV and Td ec ~ 308 
MeV, i.e. still (slightly) different transition tempera- 
tures. This gap reduces further when the corresponding 
transition temperatures from the chiral susceptibilities 
Xr/T 4 and \R are considered, see table [TT1 

In general, our results in the infinite vol- 
ume/continuum limit agree qualitatively with the 
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FIG. 4: Left diagram: Angular dependence of the quark condensate evaluated at two different quark masses and in the chiral 
limit at T = 400 MeV. Right diagram: Temperature dependence of the half-width L of the plateau in the chiral quark condensate 
in units of 7r (i.e. L/n = 1 means that the plateau extends over all angles <p = . . . 2tt) and temperature dependence of the 
chiral quark condensate A(cp) at periodic boundary conditions ip = 0. 



ones on a torus and almost quantitatively as concerns 
the transition temperatures. The most prominent 
remaining discrepancy is seen in the relative decrease 
of the ordinary condensate and the relative increase of 
the dual condensate below and above the transition 
temperature. This difference is a cutoff effect, as can 
be inferred from comparing results for the two different 
lattice spacings in the left plot of Fig. [3l In particular 
we found that cutoff effects in the quark-DSE only 
play a very minor role, whereas the cutoff effects in the 
condensate equation (JTBJ) are most significant. These 
effects will be explored further in future work. 

In general we wish to emphasize that the present cal- 
culation, although carried out with quenched lattice re- 
sults for the gluon propagator, is in itself not strictly 
quenched: our ansatz for the quark-gluon vertex is too 
simple to strictly represent the quenched theory. This 
can be seen from the fact that the dressed Polyakov-loop 
is not strictly zero below the deconfinement transition. 
Consequently we do not observe the second order decon- 
finement phase transition expected from quenched SU (2) 
Yang-Mills theory but a, more or less rapid, crossover at 
Tdec ~ 308 MeV 1 . On the other hand, as concerns the 



T 






T- dec 


301(2) 


304(f) 


305(f) 


308(2) 



TABLE II: Transition temperatures in the infinite vol- 
ume/continuum limit for the different definitions Eqs. (II3|1 . 
(|14|) and (|f 5|) . These values have been extracted from numer- 
ical results evaluated using temperature steps of one MeV 
(not shown in Fig[3}. 



chiral properties of our model, a prominent feature of 
the quenched theory not reproduced by our framework 
is the appearance of quenched chiral logarithms in the 
chiral condensate. These are well-known to be generated 
by r{ hairpin diagrams, which are not represented by our 
vertex ansatz. For the present investigation this is more 
an advantage than a drawback. Quenched chiral loga- 
rithms are most notable in the chiral limit, where they 
lead to a singularity in the chiral condensate. Since we 
do not encounter this singularity we are in a position to 
investigate both, the ordinary and the dual condensate 
in the chiral limit. 



B. The chiral limit 

In Fig. [3] we study the dependence of the quark con- 
densate on the boundary angle ip. In the left diagram 
we compare the angular dependence of the condensate 
at T = 400 MeV for two different quark masses and 
in the chiral limit. As already noted in Ref. [l(| we 
clearly see a broadening in the central dip of the graphs 
with decreasing quark mass. This can be readily under- 
stood from the loop expansion of the quark condensate, 
Eq. ([2]), which is applicable at finite quark masses and 
under presence of an ultraviolet regulator for the con- 
densate, cf. section MI B[ At sufficiently large quark 
masses large loops are suppressed by powers of 1/m. As 
a result only loops winding once around the torus should 
contribute in Eq. @ and the resulting angular behavior 
of the condensate should be proportional to cos(</?). In- 
deed, this is what we see: the result for our largest quark 



1 Note that even if our vertex were strictly quenched it is not 
clear whether the lattice input for the gluon propagator is precise 



enough to allow for an observation of the second order phase 
transition expected in SU (2) Yang-Mills theory 



9 




mass can well be fitted by only few terms in an expan- 
sion A (if) = X)^=o a ™ cos ( n f) an d the first term is by far 
the largest contribution. For smaller quark masses we ob- 
serve also sizeable contributions from terms cos(nip) with 
n > 1. In the plot, these contributions are responsible 
for the flat area around the antiperiodic boundary angle 
ip = 7r. Approaching the chiral limit this area becomes 
flatter and finally develops a derivative discontinuity at 
two finite values of if = tt ± L. These indicate the break- 
down of the loop expansion Eq. (|2|) in the chiral limit. 
Note that this is also the limit where the condensate is 
free of quadratic ultraviolet divergencies and the contin- 
uum limit can be taken; another reason why the loop 
expansion Eq. ^ is no longer applicable. 

In the plot on the right hand side of Fig. 0] we show 
the temperature dependence of the half-width L of the 
plateau in the chiral quark condensate in units of n. The 
width of the plateau is zero below the critical tempera- 
ture T c and rises quickly above, spreading up to L/ir = 
0.6 for temperatures in the region T ps 2T c . This mono- 
tonic rise then slows down considerably and seems to con- 
verge to a finite value smaller than one. However, from 
the presently available results up to T k 2.5T c one can- 
not completely exclude that the plateau finally extends 
over the whole range of boundary angles if — . . . 2-7T for 
T — > oo. 

The behavior of the chiral quark condensate at peri- 
odic boundary conditions, A(if = 0), is also a monotonic 
function of temperature as shown in the same plot. From 
a naive dimensional analysis of the condensate one may 
expect that it rises proportional to the third power of the 
temperature. This is however not correct as explicitly 
shown in our analytical analysis in appendix |SJ There 
we demonstrate that 

A V=0 (T)~T 2 for T>T C . (17) 

To extract this behavior also from our numerical data 



we fitted a polynomial of degree three in T: the result, 
A V=0 (T > T c ) = 10- 5 + 0.045(l)T+0.73(l)T 2 + lCr 7 T 3 , 
for temperatures T > 2T C clearly agrees with the ana- 
lytical analysis. This scaling behavior also agrees with 
corresponding results of Ref. fl9j . In addition it may be 
interesting to note that this T 2 -scaling of the condensate 
at periodic boundary conditions agrees with the naive 
scaling of the chiral condensate of the resulting three- 
dimensional theory in the infinite temperature limit. 

In the diagram of the left hand side of Fig. [5] we plot 
the temperature dependence of the ordinary and the chi- 
ral condensate as well as their temperature derivatives 
in the chiral limit. One clearly observes a second order 
chiral phase transition. The temperature derivative of 
the ordinary condensate diverges at the corresponding 
critical temperature T c = r = T XR / T 4 = 298(1) MeV. 
Above T c the ordinary condensate is strictly zero. For 
the deconfinement transition we observe a qualitatively 
similar behavior as for the case with finite quark mass 
discussed in the last subsection. The corresponding tran- 
sition temperature, however, has moved to a lower value, 
Tdec = 299(3) MeV, and coincides almost precisely with 
the chiral critical temperature. This is in agreement with 
corresponding extrapolations within lattice QCD [32| . 
We regard this as one of the main results of the present 
investigation. 

A combined plot of the angular and temperature de- 
pendence of the chiral limit condensate is shown in the 
right hand diagram of Fig. O We clearly observe a differ- 
ent evolution of the quark condensate at different bound- 
ary angles if. Whereas at periodic boundary conditions 
ip = 0, 2-7T the condensate is monotonically rising, cf. the 
curve in the right diagram in Fig. 0] it clearly shows sig- 
nals for a phase transition for boundary angles close to 
antiperiodic boundary conditions if = n, cf. the corre- 
sponding curve in the left diagram of Fig. Indeed for 
angles around if — tt one could (formally) define a if- 
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FIG. 6: Left diagram: A 3d-plot of the angular and temperature dependence of the quark scalar dressing function evaluated 
at zero momentum and lowest Matsubara frequency. Note that for T > 400 MeV (not shown in the plot) the scalar quark 
dressing function at periodic boundary conditions starts to rise linearly with temperature, B v =o :P =o{T) ~ T. Right diagram: 
Comparison of the temperature dependence of the corresponding dual scalar quark dressing and the dual condensate. 



dependent transition temperature that is related to the 
width L of the plateau. 

Finally we would like to make two comments. First, 
the corresponding dependence of the massive quark con- 
densate as function of (T, ip) is similar to that of its 
chiral counterpart. The most pronounced difference is 
that there is no plateau around ip = 7r, but instead a 
cosine-shaped dip as can be inferred from the left di- 
agram in Fig. 2J Second, we also determined the chi- 
ral condensate using the simplified, temperature inde- 
pendent quark-gluon interaction, Eq. ([12]). and a gluon 
propagator fixed at T = 0. This form of the interac- 
tion allows for temperature effects in the quark propaga- 
tor only, while retaining the chiral critical temperature 
t w 300 MeV. The corresponding result is qualitatively 
similar to the one with full temperature dependence in 
the interaction. However, both transitions are much less 
pronounced. The chiral condensate starts to decrease and 
the dual condensate starts to increase at much lower tem- 
peratures. Correspondingly we obtain broad curves for 
the temperature derivatives. We conclude that a quan- 
titative study of the details of the phase transitions of 
QCD has to take into account the temperature effects in 
the interaction. 



C. The dual scalar quark dressing function 

We now come back to the dual scalar quark dressing 
introduced by Eq. Our numerical results for the an- 
gular and temperature dependence of the scalar quark 
dressing function B(Q) is shown in the left diagram of 
Fig. [5J The resulting pattern is similar to the one for 
the quark condensate, cf. Fig. in particular as con- 
cerns the opening of a zero plateau around antipcriodic 



boundary conditions above the critical temperature. For 
temperatures below T c we observe a more distinct ip de- 
pendence as for the dual condensate, again indicating the 
non-quenched nature of our model for the quark-gluon 
vertex as discussed above. Above T c the scalar quark 
dressing function £^=0(0) at periodic boundary condi- 
tions stays constant for a while but then starts to rise for 
temperatures above 2T C (not shown in the 3d-plot). Our 
fit to the data leaves us with 

Bp = o lP =o(T) ~ VT for T > T c (18) 

in agreement with our analytical analysis in appendix [Bl 
and corresponding results of Ref. [r| . 

The data for the dual scalar quark dressing are shown 
on the right hand diagram of Fig. [5] and compared to the 
dual condensate. Again, the temperature dependence of 
both quantities is very similar. The temperature deriva- 
tives clearly peak at precisely the same transition tem- 
perature, indicating the deconfmement transition. Also, 
the rate of change at and around the transition tem- 
perature is quantitatively similar, although the signal 
is slightly stronger for the dual scalar dressing function 
S b ■ From a systematic perspective the dual condensate, 
or dressed Polyakov loop, has the advantage of a direct 
connection to the ordinary Polyakov loop in the limit of 
static quarks. However, the dual scalar quark dressing 
function has the advantage of being a well defined quan- 
tity in the continuum limit. After renormalization of the 
Dyson-Schwinger equation of the quark propagator this 
quantity stays finite when the cutoff is sent to infinity 
also for finite quark masses. Both quantities are promis- 
ing candidates for a further study of the deconfmement 
transitio of QCD. 
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V. SUMMARY AND CONCLUSIONS 

In this work we have investigated the chiral and the 
deconfinement transition of quenched QCD. From the 
temperature dependent quark propagator we extracted 
the order parameter for the chiral transition, the quark 
condensate, and order parameters for the deconfinement 
transition, the dressed Polyakov loop and the dual scalar 
quark dressing. We determined the quark propagator 
from its Landau gauge Dyson-Schwinger equation in the 
imaginary time Matsubara formalism. We emphasize 
once more that our truncation scheme, introduced in 
Ref. [l6| , implements a realistic temperature dependence 
of the gluon propagator and takes into account temper- 
ature effects in the quark-gluon vertex. It therefore goes 
beyond the simple models used in previous works. 

When comparing results from the infinite vol- 
ume/continuum DSEs with corresponding ones on a 
torus we found interesting effects due to different renor- 
malization prescriptions for the quark propagator: the 
deconfinement transition temperature on a torus is 
shifted to smaller temperatures, i.e. from Tdec ~ 320 
MeV to Tdec ~ 308 MeV when artificial temperature ef- 
fects in the renormalized quark mass are avoided. This 
value agrees precisely with Td ec ~ 308(2) MeV in the in- 
finite volume/continuum limit reflecting the fact that the 
remaining cutoff artefacts on our torus do not affect the 
transition temperature. The same is true for the chiral 
transition temperature. We find r = 301(2) MeV from 
the temperature derivative of the chiral condensate and 
T Xb jt 4 = 304(1) MeV from the normalized chiral sus- 
ceptibility. In summary we conclude that the chiral and 
deconfinement transition temperatures are only slightly 
different for finite quark masses. This finding will be fur- 
ther tested in future work using more precise data for the 
gluon propagator. 

In the chiral limit this agreement becomes exact 
(within numerical errors). We find a second order chiral 
phase transition at T c = 298(1) MeV and a similar tem- 
perature for the deconfinement transition, Td ec = 299(3) 
MeV. It is worth to emphasize again that both transi- 
tion temperatures are extracted from the properties of 
the quark propagator and therefore related to the under- 
lying properties of the Dirac operator. Therefore in a 
sense discussed in more detail in Refs. @, H, [1(1 [ll|, Ojj] 
the chiral and the deconfinement transition arc closely 



connected. 

The framework used in this work is quenched SU (2) 
Yang-Mills theory. Our transition temperatures may 
be translated into the corresponding ones of quenched 
SU{3) QCD using the relations T c /^ = 0.709 (SU(2)) 
and T c /y/a = 0.646 (SU(3)) between the respective crit- 
ical temperatures and the string tension [331 ] . The re- 
sulting transition temperatures are then T XR /T i ~ 277 
MeV and Td ec ~ 281 MeV for the massive case and 
T Xr/T 4 ss T dec ps 272 MeV in the chiral limit. In or- 
der to work in the full, unquenched theory we would 
have to take into account quark-loop effects in the gluon 
propagator and meson effects in the quark-gluon vertex 
[27j . These effects will shift the transitions tempera- 
tures below T = 200 MeV, see 0, H] for latest results 
for Nf = 2 + 1 quark flavors. As concerns the dual con- 
densate and scalar dressing function in the unquenched 
formulation one needs to carefully take into account ef- 
fects due to the Roberge- Weiss symmetry [2(J. This is 
because of the formal similarity of the continuous bound- 
ary conditions for the quark field to an imaginary chem- 
ical potential, see [l!| for details. An investigation of 
these effects as well as a detailed study of volume and 
discretization artefacts is deferred to future work. 
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APPENDIX A: EXPLICIT FORM OF THE DSES 
FOR THE QUARK DRESSING FUNCTIONS 

DSEs in the infinite volume/continuum limit: 

The explicit form of the DSEs in the infinite vol- 
ume/continuum limit in the truncation specified in sec- 
tion IIII Al is given by 
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UJ 2 C 2 (q)+cpA 2 (q) + B 2 (q) 
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2 A 
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» T A (q, fc) + A L (p) (r A (q, fc) -f - r c (g, fc) ^) 
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(A3) 



with q = (g, u> q ), k = (fc, Wfc), p = (p, cj p ) = g — fc, and the gluon propagator 



A T (p) 



-Mp) 



A l (p) 



(A4) 



and the vertex dressing functions 



r°(g,fc)=r(g-fc) 



C(q) + C(k) 



and 



r A (g,fc)=r(g-fc) 



A(q)+A(k) 



(A5) 



with 



r(g) 



dx + g 2 



d 2 + g 2 A 2 + g 2 



/3 a(Ai)ln[g 2 /A 2 + l] 



-in 



(A6) 



r 



DSEs on a torus: 

The corresponding form of the DSEs on a torus is simi- 
lar to the one given in Eq. (|ATj) ,Eq. ljA"2]) and Eq. (|A"3"]) 
with the exception that the three-dimensional momen- 
tum integral is replaced by a sum over momenta. In the 
three Cartesian spatial directions we thus have Matsub- 
ara sums counting momenta p n = J2i=i / L){rii)ei, 
where e; are Cartesian unit vectors in Euclidean momen- 



tum space. In 31| a technique of rearranging the three 
Cartesian sums into a sum over hypcr-sphcres and a sec- 



ond sum over the individual momenta on each hyper- 
sphere has been developed that is convenient for the nu- 
merical treatment of the DSEs. This sum is given by 



L 3 



(A7) 



where j counts spheres with p 2 = const, and m numbers 
the grid points on a given sphere. The corresponding 
momentum vectors are denoted by p^ m - 



13 



APPENDIX B: LARGE TEMPERATURE 
SCALING OF B v=0 AND <V>V>) v =o 



This result is again in excellent agreement with our nu- 
merical results discussed in section HVBI 



In this section we analyze the large temperature scal- 
ing of the scalar quark dressing function B v= o(p, T) and 
the quark condensate (ipip) v= q(T) with temperature. To 
this end we first analyze the large temperature limit of 
the DSE Eq. (|A2[) for the scalar quark dressing func- 
tion. With periodic boundary conditions ip = only the 
zcroth Matsubara frequency uj q = contributes in the 
Matsubara sum; all others are suppressed by powers of 
the temperature T. Furthermore at large enough tem- 
peratures the temperature effects in the gluon dressing 
functions and the vector dressing functions A and C of 
the quark propagator and quark-gluon vertex are arbi- 
trary weak and therefore can be neglected. Working in 
the chiral limit we then have 



B(k,T) ~ T 



T / dq 



d 3 q B(q, T) 3T{q-k) 
B 2 (q,T) (q-_jfe)2 

B(q,T) 



(2tt) 3 q 2 



(Bl) 



B*{q,T) ' 

The second line follows from the first line because of the 
temperature independence of the terms in the angular 
integral. The resulting integral in Eq. (|B1[) is dominated 
from low momenta q 2 < B 2 (q, T) and consequently the 
integrand scales like l/B(q,T) leading to 



B(k,T) ~ VT 



for T>T C; 



(B2) 



in agreement with our numerical results discussed in sec- 
tion [IVC] and corresponding results in Ref. [lji[. We also 
confirmed that indeed the scaling of B(k,T) with T is 
independent of the three-momentum k. 

Having identified the large temperature scaling of 
B(k,T) we now analyze the integral for the quark con- 
densate Eq. (flU)) given by 

m d3q 



B 



T / dqq 



-a 



(2tt) 3 u 2 C 2 + q^A 2 + B 2 
B(q,T) 



B 2 (q,T) 



(B3) 



Here the second line is obtained by again neglecting the 
temperature scaling of the vector dressing functions A, C 
and the contributions of all but the lowest Matsubara 
frequency. The resulting integral has an extra factor of 
q 2 as compared to the one in Eq. (|B 1[) . It is therefore 
not dominated by momenta smaller than B 2 . Recalling 
that B(0, T) ~ B(q, T) ~ \Ff we proceed by the integral 
transformation q — > q'B(Q,T) which leads to 



{il>il>) v =o 



T I dq'q' 2 B 3 (Q,T) x 



B(q,T) 



B 2 {q,T) (e 2 B 2 (0,T)/B 2 (q,T) + \ j 
T 2 , for T>T C . (B4) 



APPENDIX C: SCALING VS. DECOUPLING 

In this appendix we demonstrate that the deep infrared 
behavior of the gluon propagator and the quark-gluon 
vertex is not related to the details of the deconfincment 
transition at Tdec- The reason why we are interested 
in this question is an ongoing discussion over the past 
years on the status of two possible types of analytical 
infrared solutions of Yang-Mills theory. These two types 
of solutions arc characterized by a one parameter family 
of boundary conditions on the dressing function of the 
ghost propagator, G(p 2 ), at zero momentum. If G(p 2 = 
0) = oo one observes the 'scaling' solution of infrared 
Yang-Mills theory, characterized by power laws in all one- 
particle irreducible Green's functions of the theory [HI, 
[H, [HI]. The other possibility, G(p 2 = 0) < oo, describes 
'decoupling' characterized by a finite gluon propagator 
and a finite ghost dressing function at zero momenta. [HI, 
H3, HH- Only the first of these two possibilities, scaling, 
is in agreement with the Kugo-Ouma criterion of well- 
defined global color charges [39|, |40| ■ 

Nevertheless, both possibilities satisfy a confinement 
criterion developed in Ref. .lj§]: both types of infrared 
solutions lead to a confining Polyakov-loop potential. It 
has also been noted in Ref. [l8|] that the confinement- 
deconfinement phase transition extracted from the tem- 
perature dependence of the Polyakov-loop potential is 
not qualitatively affected by the choice of scaling or de- 
coupling. Here we demonstrate that the same is true for 
the deconfincment transition from the dual quark con- 
densate. 

To this end we note that the choice of the fit-function 
Eq. (fTX))) for the gluon propagator and the ansatz for the 
quark-gluon vertex Eq. (|11[) correspond to 'decoupling' 
in the deep infrared. An alternative choice that fits the 
lattice results of Ref. [I3| for the temperature dependent 
gluon propagator equally well is given by 



7 scaling / -> 



(q,w q ,T) = Z T {q,uJ q ,T) 



Z L (q^ q ,T) = Z L (q,uj q ,T), 



2K+1/2-1 

<" + A! (C1) 



(C2) 



where Zr,L{q, w q ,T) are given by Eq. (fTX)]) . The effect 
of the extra factor is to convert the infrared behavior 
ZT(q,uj q ,T) ~ q 2 of the decoupling type of gluon dress- 
ing function to the scaling behavior Z^ ahn9 (q,oj q ,T) ~ 
(q 2 ) 2K+1 / 2 with the anomalous dimension n w 0.40 [4l| of 
the three-dimensional theory. Corresponding scaling so- 
lutions from DSEs have been reported in Ref. JlT\ . Note 
that the electric part of the gluon propagator always de- 
velops a screening mass and is therefore not modified in 
Eq. (|C2[) . The scale Ajr is typically much smaller than 
Aqcd', here we choose A IR = 100 McV in agreement 
with the results of 23]. The scaling of the quark-gluon 
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FIG. 7: Temperature dependence of the dressed Polyakov-loop Ei and the conventional quark condensate = {tptp)^^ 
together with their derivatives for m = lOMeV in the infinite volume/continuum limit. The thick lines, already shown in 
the right diagram of Fig. [3] correspond to the decoupling type of interaction whereas the thin lines are the scaling type of 
interaction Eqs. (fClj) , (|C2) and (JClJ. 



vertex in the quenched theory is also given by a power 
law [12] and can be represented by the ansatz 

C"" 9 (g^,p) = r v ( q ,k,p) ( g2 f A ,J ' (C3) 

where q is the gluon momentum and T v (q, k,p) is given 
by Eq. (HH). 

Our numerical solutions for the quark condensate and 
the dressed Polyakov loop are given in Fig. [71 The re- 
sults with decoupling type of interaction already given in 
Fig. [3] in the main body of this work are shown together 
with the results obtained from the scaling type of inter- 



action. The differences are tiny. Below the transitions 
the effects are not visible. Here the dynamically gener- 
ated quark mass in the quark-DSE shields the interaction 
in the deep infrared; the results are independent of the 
behavior of the interaction below A//j. We also find the 
same transition temperatures for both types of solutions. 
In the high temperature phase where chiral symmetry is 
partly restored we observe small effects in the ordinary 
condensate and almost no effect in the dual chiral con- 
densate. The main aspects of the two transitions are not 
concerned. Whether there is some impact on the proper- 
ties of quarks within the quark-gluon plasma remains to 
be investigated. 



[2] 
[3] 

[4] 

[5] 

[6] 
[7] 

[«] 
[9] 
[10] 

[11] 



A. M. Polyakov, Phys. Lett. B 72 (1978) 477; 

L. Susskind, Phys. Rev. D 20, 2610 (1979). 

F. Karsch and E. Laermann, arXiv:0305025 [hep-lat]. 

Y. Aoki, Z. Fodor, S. D. Katz and K. K. Sza bo, Phys. 

Lett. B 643 (2006) 46 [arXiv: hep-lat /0609068] . 

A. Bazavov et al, Phys. Rev. D 80 (2009) 014504 

[arXiv:0903.4379l [hep-lat]] . 

Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, 
S. Krieg and K . K. Szabo, JHEP 0906, 088 (2009) 
[arXiv:0903.4T55l [hep-lat]] . 

J. Greensite, Prog. Part. Nucl. Phys. 51 (2003) 1 



arXiv:hep-lat/0301023| . 

Gattnar, C. Gattringer, K. Langfeld, H. Reinhardt, 
A. Schafer, S. Solbrig and T. Tok Nucl. Phys. B 716, 
105 (2005) [arXiv:hep-lat/0412032] . 
C. Gattringer, Phys. Rev. Lett. 97 (2006) 032003, 
[arXivihep-iat /0605018"] . 

F. Bruckmann, C. Gattringer and C. Hagen, Phys. Lett. 

B 647 (2007) 56, [arXiv:hep-lat/0612020| . 

F. Synatschke, A. Wipf and C. Wozar, Phys. Rev. D 75 

(2007) 114003, [arXiv:hep-lat/0703018l . 

F. Synatschke, A. Wipf and K. Langfeld, Phys. Rev. D 

77 (2008) 114018, [arXiv:0803.0271l [hep-lat]].' 



[12] E. Bilgici, F. Bruckmann, C. Gattringer and C. Hagen, 

Phys. Rev. D 77 (2008) 094007, [arXiv:0801.405T1 . 
[13] F. Bruckmann, C. Hagen, E. Bilgici and C. Gattringer, 

PoS LATTICE2008 (2008) 262, |arXiv:0810.0899] . 
[14] J. Danzer, C. Gattringer and A. Maas, JHEP 0901 

(2009) 024, i arXiv:0810.3973l [hep-lat]] . 
[15] E. Bilgici, F. Bruckmann, J. Danzer, C. Gattringer, 

C. Hagen, E. M. Ilgenfritz and A. Maas. larXTv:0906.3957l 

[hep-lat] . 

[16] C. S. Fischer, P hys. Rev. Lett. 103, 052003 (2009), 

larXiv:0904.2700l [hep-ph] . 
[17] A. Cucchieri, A Maas and T. Mendes, Phys. Rev. D 75, 

076003 (2007), [arXiv:hep-lat/0702022l . 

[18] J. Braun, H. Gies and J. M. Pawlowski. larXiv:0708.24l3l 
[hep-th]; F. Marhauser and J. M. Pawlowski, 
larXiv:0812.lT44l [hep-ph] . 
[19] J. Braun, L. Haas, F. Marhauser and J. M. Pawlowski, 

larXiv:0908.0008l [hep-ph] . 
[20] A. Roberge and N. Weiss, Nucl. Phys. B 275 (1986) 734. 
[21] R. Williams, C. S. Fischer a nd M. R. Pennington , 
Phys. Lett. B 645 (2007) 167 [arXiv:hep-ph/06120611 

[22] 



arXiv:0704.2296 [hep-ph] 
C. D. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. 



15 



45, SI (2000), [arXiv:nucl-t h/0005064]. 

[23] C. S. Fischer, A. Maas and J. M. Pawlowski, Annals 
Phys. 324, 2408 (2009) [arXiv:0810.1987l [hep-ph]]. 

[24] A. Maas, J. Wambach and R. Alkofe r, Eur. Phys. J. C 
42 (2005) 93, [arXiv:hep-ph/0504019| ; 

[25] A. Sternbeck, E. M. Ilgenfritz, M. Muller-Preussker, 
A. S chiller and I. L. Bogolu bsky, PoS LAT2006 (2006) 
076 [arXiv:hep-lat/0610053| . 

[26] J. S. Ball and T. W. Chiu, Phys. Rev. D 22 (1980) 2542. 

[27] C. S. Fi scher and R. Wil liams, Phys. Rev. D 78 (2008) 
074006, |arXiv:0808.3"372l [hep-ph]] . 

[28] A. Bender, D. Blaschke, Y. Kalinovsky and 
C. D. Roberts, Phys . Rev. Lett. 77 (1996) 3724, 
[arXiv:nucl-th/9606006] . 

[29] P Maris, C. D. Roberts, S. M. S chmidt and P. C. Tandy , 
Phys. Rev. C 63, 025202 (2001) larXiv:nucl-th/0001064| . 

[30] D. Horvatic, D. Klabucar and A. E. Radzhabov, Phys. 
Rev. D 76 (2007) 096009 [arXiv:0708. 12601 [hep-ph]]. 

[31] C. S. Fischer, R. Alk ofer and H. Reinhardt, Phys. Rev. 
D 65 (2002) 094008, [arXiv:hep-ph/0202195l ; C. S. Fis- 
cher an d M. R. Pennington, P hys. Rev. D 73 (2006) 
034029, [arXiv:hep-ph/0512233| ; C. S. Fischer, A. Maas, 
J. M. Pawlowski and L. von Smekal, Annals Phys. 322 
(2007) 291 6. [arXiv:hep-ph/070T050] . 

[32] F. Karsch, [arXiv:hep- lat /9903031 1 

[33] J. Fingberg, U. M. Heller and F. Karsch, N ucl. Phys. 
B 392 (1993) 493; [arXiv:hep-lat /9208012 j . O. Kacz- 



[34] 
[35] 

[36] 
[37] 

[38] 
[39] 

[40] 



Lett. B 543 (2002) 41; [arXiv:hep-lat 70207002] . B. Lu- 
cini, M. Teper and U. Wenger, JHEP 0502 (2005) 033; 
R. Alkofer, C. S. Fischer and F. J. Llanes-Estrada, Phys 
Lett. B 611, 279 (2005) [arXiv:hep-th/0412330| . 
C. S. Fischer and J. M. Pawlowski. Phys. Rev 



Pawlowski, Phys. Rev. D 75, 
025012 (2007) larXiv:hep-th /0609009] ; Phys. Rev. D 80, 
025023 (2009) |arXiv:0903 21931 [hep-th]]. 
R. Alkofer, M. Q. Huber and K. Schwenzer, 
larXiv:0801.2762l [hep-th] ; 

Ph. Boucaud, J. P. Leroy, A. L. Yaouanc, J. Micheli, 
O. Pen e and J. Rodrigu ez-Quintero, JHEP 0806, 012 
(2008) |arXiv:0801.272"Tl [hep-ph]]. 

A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. 
D 78, 025010 (2008) |arXiv:0802. 18701 [hep-ph]]. 
T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66, 1 
(1979) [Erratum Prog. Theor Phys. 71, 1121 (1984)]; T. 
Kugo, ]arXiv:hep-th/9511033| 

C. Lerche and L. von Smekal, Phys. Rev. D 65, 
125006 (2002) [arXiv:hep-ph/0202194l ; J. M. Pawlowski, 

D. F. Litim, S. Nedelko and L. von Smekal, Phys. Rev. 
Lett. 93, 152002 (2004) [arXiv:hep-th/0312324] . 

(2002) 



[41] D. Zwanziger, Phys. Rev. D 65, 094039 
[42] 



[arXiv:hep-th/0109224| ; 

R! Alkofer, U. S. Fischer, F. J. Llanes-Estrada 
and K. Sc hwenzer, Annals Phys. 324, 106 (2009) 
[arXiv:0804.3042l [hep-ph]] ; 



marek, F. Karsch, P. Petreczky and F. Zantow, Phys. 



